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ABSTRACT 

We report on our second-year campaign of X-ray follow-up observations of uniden- 
tified Fermi-LAT 7-ray sources at high Galactic latitudes > 10°) using the X-ray 
Imaging Spectrometer onboard the Suzaku X-ray Observatory. In this second year of the 
project, seven new targets were selected from the First Fermi-LAT Catalog, and studied 
with 20 — 40ks effective Suzaku exposures. We detected an X-ray point source coincident 
with the position of the recently discovered millisecond pulsar PSR J2302+4442 within 
the 95% confidence error circle of IFGL J2302. 8+4443. The X-ray spectrum of the de- 
tected counterpart was well fit by a blackbody model with temperature of kT ~ 0.3 keV, 
consistent with an origin of the observed X-ray photons from the surface of a rotating 
magnetized neutron star. For four other targets which were also recently identified with 
a normal pulsar (IFGL J0106. 7+4853) and milhsecond pulsars (IFGL J1312. 6+0048, 
J1902.0-5110, and J2043. 2+1709), only upper limits in the 0.5 - lOkeV band were ob- 
tained at the flux levels of ~ 10~^^ ergcm"^ s~^. A weak X-ray source was found in the 
field of IFGL J1739. 4+8717, but its association with the variable 7-ray emitter could 
not be confirmed with the available Suzaku data alone. For the remaining Fermi-LAT 
object IFGL J1743.8— 7620 no X-ray source was detected within the LAT 95% error 
ellipse. We briefly discuss the general properties of the observed high Galactic-latitude 
Fermi-LAT objects by comparing their multiwavelength properties with those of known 
blazars and millisecond pulsars. 
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8 Subject headings: Stars: pulsars: general — Stars: pulsars: individual (PSR J2302+4442) — 

9 Galaxies: active — Gamma rays: general — X-rays: general 



1. Introduction 



20 



23 



Since its success ful launch in 2008 J une, the Large Area Telescope (LAT) onboard the Fermi Gamma- 
ray Space Telescope (lAtwood et al.ll2009) has enabled many important breakthroughs in the understanding 
of the origin of high energy 7-ray emissions of various classes of astrophysical objects. The number of 
detected 7-ray sourc es increased dramatically, from 271 objects listed in the 3rd EGRET C atalog (3EG; 



Hartmanetal] 1 199911^ to 1873 in the Second Fermi-LAT Catalo g (2FGL; 



Abdo et al 



20111). About 800 



7-ray sources included in 2FGL have been identified as blazars jAckermann et al. 2011 ). i.e., jetted ac 



17 tive galactic nuclei (AG N) characterized by strong relati yistic beaming. Ot her associations inclu ded pulsars 



(e.g.jAbdo et al. I boiOal) . high-ma ss X-ray binaries (e .g., lAbdo et al. Il2009bh . radi o galaxies (e.g., 



- „ _ jAbdo et al 

2010gl) . pulsar w ind nebulae (e.g., Abdo et aPboiOe ). supe rnova remnants (e.g ., Abdo et al. 2010d[) . globu 



lar clusters (e.g., Abdo et aPboiOi ). starburs t galaxies (e.g., Abdo et al. 2010c ). and distinct objects like the 
Large Magellanic Cloud (|Abdo et al I boiohh . However, no obvious counterparts at longer wavelengths have 
been found for as much as 38% of Fermi-LAI objects so that several hundreds of GeV sources currently 
remain unassociated with any known astrophysical systems. Fortunately, an improved localization error for 

24 the Fermi-L AT (typical 95% confidence radii rgs ^ 0°.l - 0°.2, and even 0°.005 - 0°.01 for the bright- 

25 est sources; lAbdo et al.ll201ll ). when compared to that of EGRET (typical rgs ~ 0°.4 — 0°.7), allows for 



26 much more effective follow-up studies at radio, optical, and X-ray frequencies, which can help to unravel 

27 the nature of the unidentified 7-ray emitters. 

28 In this context. X-ray follow-up observations of unidentified Fermi-hAI objects are of particular im- 
portance, since some classes of astrophysical sources of 7-rays such as AGN are strong X-ray emitters 
as well, while the others like most of 7-ray emitting pulsars are faint X-ray sources. Note that assum- 
ing the keV-to-GeV emission continuum in the form of a broad-band power-law {Fi^ oc z^^""^), which 
could be a relatively good zero-order approximation in the case of blazar sources but not necessarily 
in the case of other classes of 7-ray emitters, the monochromatic X-ray flux energy density scales as 
[i^i^r.]ikev = (IkeV/O.lGeV)^"""^ x [iyF^]o,iGcV =^ 3 x lO^^ x [z^F,,]o.igcV for a relatively flat spectral 
index of a^^ ~ 0.5. Hence, if an X-ray counterpart of a bright Fermi-LAI source is characterized by, e.g., 
[i^Fj^Jo.iGeV — 10~^^ ergcm"^ s~^ and the X-ray-to-7-ray power-law emission continuum with the slope 
a^-y > 0.5, such source can be expected to be detectable with modern X-ray instruments such as Chandra, 
XMM-Newton, Swift, and Suzaku within reasonable exposure times. In particular, a point source search to 

39 the level of \u Fl,]^\(pM ^ fl()~-^'^ — 10^^^) ergcm~^ s~^ is easily attaina ble with the X-ray Imaging Spec- 

40 trometer (XIS: lKoyama et al.ll2007r ) onboard Suzaku (iMitsuda et al. l2007h with relatively short exposures o: 



See also lCasandiian & Grenien ( 120081) for the revised catalog of EGRET 7-ray sources. 
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41 few tens of ksec (e.g.. lAkamatsu et al.ll2011r) . In the case of a positive detection, correlated flux changes at 

42 X-ray and 7-ray frequencies provide an identification. The lack of correlated variability, or non-detection of 

43 an X-ray counterpart, provide on the other hand only circumstantial evidence regarding the nature of a stud- 

44 led target. Yet in many cases such evidence may be crucial, since the non-detection of an X-ray counterpart 

45 despite a long, dedicated observation of a bright Fermi-LAI object may disprove a potential association with 

46 given classes of astrophysical sources. That is because, as mentioned above, only a few established high- 

47 energy emitters are that bright in 7-ravs b ut very faint in X-rays (e.g., Geminga pulsar; see the discussion in 



48 



Thompsonll2004l : iMatsumoto et al.ll2007l) . 



49 Thus motivated we started a project to investigate the nature of unidentified high Galactic-latitude 

50 Fermi-LAI objects through deep X-ray follow-up observations with Suzaku XIS. The re sults of the first-year 

51 campaign conducted over the span of Suzaku A04 were presented in lMaeda et al.l (|201lh . The A04 program 

52 include d four steady/weak ly variable Fermi-LKI sources from the initial Fermi-LKI Bright Source List 

53 (OFGL: lAbdo et al.ll2009ah and can be summarized as follows. The X-ray counterpart for one of the bright- 

54 est unassociated Fermi-LAT objects, IFGL J1231. 1-1410 (also detected by EGRET as 3EG J1234-1318 

55 and EGR J1231 — 1412), was found. The X-ray spectrum of the counterpart was well fit by a blackbody 

56 model with a temperature of kT ~ 0.16 keV plus an additional power-law component dominating above 

57 2 keV photon energies. This power-law component was confirmed in subsequent Swift and XMM-Newton 

58 exposures. Consid ering a recent identi fica tion of IFGL J1231. 1 — 1410 with the millisecond pulsar (MSP) 

59 PSR J123 1 — 141 1 (IRansom et al.ll201lh . in lMaeda et al.l we concluded that the detected thermal X-ray pho- 

60 tons originate from the surface of a rotating magnetized neutron star, while the non-thermal X-ray compo- 

61 nent is most likely produced within the pulsar magnetosphere. In the case of IFGL J1311.7— 3429, two pos- 

62 sibly associated X-ray point sources were discovered, one of which is now excluded from the smaller error 

63 ellipse of the GeV emitter as catalogued in the 2FGL (Abdo et allboil ). The identification of the remaining 

64 X-ray counterparts with the respective 7-ray objects remain uncertain despite a robust determination of the 

65 the spectral and variability properties of the X-ray sources. In the case of IFGL J1333.2-I-5056, we found 

66 several weak X-ray sources within the 7^2rm/-LAT error circle, and speculated on the AGN nature of the 

67 target. Finally, one X-ray point source was detected at the edge of the error ellipse of IFGL J2017. 3-1-0603. 

68 The physical connection was however viewed as unlikely, since the X-ray source did not coincide with the 

69 location of the MSP PSR J20 17-1-0603 discovered by th e Nangay radio telesc ope which constituted a more 

70 highly probable association with the FennZ-LAT object dCognard et al.ll201lh . The MSP identification was 

71 later indeed confirmed by the detection of the pulsed emission in the Fermi-LAT data, with the same period 

72 as the radio pulsations. 



73 In this paper, we report the results of our second-year campaign, conducted over the span of Suzaku 

74 AOS (2010 April to 201 1 March) which included observations of seven ferm/-LAT sources located at high 

75 



Galactic latitudes (|6| > 10 °). The targets were selected from the First Fermi-LAT Catalog of point sources 



76 (IFGL; lAbdo et al.l l2010bh as objects unidentified at the time of writing of the Suzaku AOS proposal 

77 Since then however, foi ir of the selected ta rgets have been associated with MSPs: IFG L J1902.0— SllO 



78 with PSR J1902-S10S dCamilo et al.ll201lh . IFGL J 2043.2-F1709 with PS R J2043-F1711 dGuillemot et al 



201 ih . IFGL J2302. 8-1-4443 with PSR J2302-1-4442 dCognard et all 1201 ih . and IFGL J1312.6-1-0048 with 



-4- 



83 



PSR J1312+00 (lAbdoetalJbOllh . And mo reover, IFGL J0106 .7+4853 have very recendy been associ- 
ated with a normal pulsar PSR J0106+4855 dpietsch et al. Iboill) . The XMM-Newton and Sw ift satellites 
detected the weak X-ray counterpart of the MSP in IFGL J2302.8-I-4443 (Icognard eXal]|201l|). In the fol- 
lowing, our new Suzaku observations and data reduction procedure are described in section |2l The analysis 
results are given in section[3j and discussed further in section|4]in the context of multiwavelength studies of 
unidentified Fermi-LKT objects. 



87 



2. Suzaku Observations and Data Analysis 
2.1. Observations and Data Reduction 
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We observed seven unidentified Fermi-LAI sources with the Suzaku XIS. These seven targets were 
chosen from the IFGL catalog according to the following selection criteria: (i) no association claimed at the 
time of the submission of the Suzaku proposal, (ii) sources located more than ±10° away from the Galactic 
plane, and (iii) the detection significance in the LAT 7-ray band (> 100 MeV) exceeding 14cr in the IFGL 
catalog. The thus selected targets are listed in Table[T] together with the corresponding Suzaku observation 
logs and the IFGL 7-ray fluxes and photon indices. 



The observations were c onducted with three XIS detector and the Hard X-ray Detector (HXD: lKokubun et al 



20071 : iTakahashi et al.ll2007h . The XIS detector is composed of four CCD cameras. One of the four CCD 



cameras (XISl) is back-illuminated CCD and the others (XISO, XIS2, and XIS3) are front-illuminated 
CCDs. The operation of XIS2 ceased in 2006 November because of the contamination by a leaked charge. 
Since none of the studied sources have been detected with the HXD, below we describe the analysis of only 
the XIS data. The XIS was operated in pointing mode and the normal clocking mode, combined with the 
two editing modes 3x3 and 5 x 5 for five targets, and only one editing mode 3 x 3 for IFGL J2302. 8-1-4443 
and IFGL J13 12.6-1-0048 because of the telemetry limit. 

We conducted all the data reduction and analysis with HEAD AS software version 6.9 and the caUbration 
database (CALDB) released on 2010 July 30. First, we combined the cleaned event data of the two editing 
modes using xselect. Then we removed the data corresponding to the epoch of low-Earth elevation 
angles (less than 5°), as well periods (and 60 sec after) when the Suzaku satellite was passing through the 
South Atlantic Anomaly (SAA). Moreover, we also excluded the data obtained when the Suzaku satellite 
was passing through the l ow Cut-Off Rigi dity (COR) of below 6 GV. Finally we removed hot and flickering 
pixels using sisclean (|Day et al.lll998r) . 



2.2. Data Analysis 

X-ray images of each target were extracted from the two operating front-illuminated CCDs (XISO and 
XIS3). In the image correction procedure we applied a 'non X-ray backg round' subtraction, an exposure 
correction, and a vignetting correction (for details see iMaeda et al.ll201 ih . We then combined the images 
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113 from both CCDs and smoothed the thus obtained final maps using a Gaussian function with a = 0'.28. 

114 The resulting images are presented and discussed in the next section[3] Although all the 7-ray targets were 

115 initially selected from the IFGL catalog, in all the corre sponding figures, thick green ellipses denote the 



129 

130 
131 



116 more precise 95% position errors from the 2FGL catalog (lAbdo et al.ll201 ih as described in section lZT] 



117 For the further analysis, we selected source regions around each detected X-ray source within the 2FGL 

118 error ellipses. The radii of the source extraction regions, denoted in the figures below by thin green circles, 

119 were set as 1', unless otherwise stated. The corresponding background regions with radii of 3' were taken 

120 from the low count rate area in the same XIS chips (dashed green circles). We set the detection threshold 

121 for X-ray sources at Aa, based on the signal-to-noise ratio defined as a ratio of the excess events above a 

122 background to its standard deviation assuming a Poisson distribution. The X-ray source positions and the 

123 corresponding errors were estimated by 2D Gaussian fits, as summarized in Table|2] 

124 For the timing analysis, light curves from the front-illuminated (XISO, XIS3) and back-illuminated 

125 (XISl) CCDs were combined; the corresponding backgrounds were subtracted using Icmath. The light 

126 curves constructed in this way provide the net-count rates. To quantify possible flux variations, the test 

127 was applied to each light curve using Icstats. For the X-ray spectral analysis, we generated the RMF files 

128 for the detector resp onse and the ARF files for the effective area using xisrmf gen and xissimarf gen 
( Ishisaki et al In order to improve the statistics, we added X-ray counts from the two front-illuminated 



CCDs, using mathpha with no error propagation so that the resulting data follow a Poisson distribution, 
and then combined the response files using the marf rmf and addrmf commands. In the case of the 7-ray 

132 targets with no detected X-ray counterparts within the 2FGL error ellipses, we calculated 90% confidence 

133 level flux upper limits at the positions of the 7-ray emitters assuming an absorbed power-law model. Un- 

134 certainties of the model spectral parameters are computed at the 90% confidence levels. The results of the 

135 spectral fitting are summarized in Tabled and discussed below in more detail. 



3. Results 



In this section, we first present the results of the new Suzaku observations in the order of R.A. for the 
sources for which we detected X-ray counterparts, followed by the analysis results of the remaining tar- 
gets with no detected X-ray counterparts. All seven targets have also recently been observed by Swift, and 
below we briefly compare the results of the Suzaku and Swift observations. We search for radio, infrared 
and o ptical counterparts for all the detected X-ray sources using the NRAO VLA Sky Sur vey (NVSS) cat 



142 alog (ICondon et al 



19981) . the Two Micron All Sky Sur vey (2MASS) point source catalog (ISkrutskie et al. 



2006I) . and the USNO-Bl.O catalog jMonet et al. Iboosi) . 



3.1. IFGL J1312.6+0048 



145 

146 



We discovered one X-ray point source inside the 2FGL error ellipse of IFGL J1312. 6-1-0048, at [RA, 
Dec]=[198°.235(2), 0°.835(2)], designated hereafter as Suzaku J1312-I-0050, with a detection significance 
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148 



of 9a (total of 178 net source co unts from three detect ors). This X-ray source has recently been detected also 
by Swift X-ray Telescope (XRT: burrows et al. Iboosb at [RA, Decl^ [265.4229, 0.8348] with 90% position 



error ra dius rQn% = 6".S. Since the position accuracy of Swift XKm 5") is better than that of Suzaku XIS 



150 (~ 19^^ lUchiyama et al. II2008I) . we searched for optical and radio counterparts of the X-ray source using 

151 the Swift position. We found one possibly related optical emitter USNOB 0908-0218088 at about 5"away 

152 from the X-ray source ([RA, Dec] = [198.241972(7), 0.83448(5)], and magnitudes B2 = 19.80, R2 = 19.32, 

153 and I = 18.11), while no radio and infrared counterparts were discovered. 

154 The light curve of Suzaku J1312-I-0050 with a time bin of 5774 s and its spectrum are presented in 

155 Figures|2]and|3l respectively. The X-ray spectrum when fitted with an absorbed power-law model returned 

156 negligible hydrogen column density. We therefore fixed iVn = and repeated a power-law fit obtaining 

157 the photon index of F = 1.9^q 3 with x^/d.o.f.=9.9/13. The derived energy flux in the 2 — lOkeV photon 

158 energy range is 8.OI26 ^ 10~^^ ergcm"^ s~^. To test for flux variability, we performed a fit with a 

159 constant count rate of 1.02 x 10^^ ct s^^ resulting in x^/d.o.f. = 14.5/7. This indicates that the X-ray source 

160 is variable on the timescale of a few hours, with a probability of ~ 96%. 



The targeted 7-ray object has recently been associated with the MSP PSR J13 12-1-00 (lAbdo et al.ll2011h . 
Unfortunately, the position of the radio pulsar is still not available publicly, and hence at this moment we 
cannot claim nor reject the coincidence of PSR Jl 3 12-1-00 with the detected X-ray source. The variability 
revealed by the Suzaku data, along with the presence of an optical counterpart, implies however that the de- 
tected X-ray source is not likely associated with the MSP PSR J 13 12-1-00. For this reason we have evaluated 
the 90% confidence upper limit to the 0.5 — 2keV and 2 — lOkeV fluxes at a random position within the 
2FGL error ellipse of IFGL J1312.6-I-0048 (excluding Suzaku J1312-I-0050), which should correspond (at 
least roughly) to the upper limits for the X-ray emission of PSR J13 12-1-00. Assu ming an absorbed po wer- 
law model with nn = 2.1 (this value is taken from the LAB Survey of Galactic HI: iKalberla et al. 1 12005b and 



photon index F = 2, these read as < 1.3 X 10 ^^ergcm ^s ^ (0.5 — 2keV) and < 9.3x 10 ^^ergcm ^s 
(2 - lOkeV), respectively. The imphed ratio of the 0.1 - 100 GeV and 2 - lOkeV flu xes F^/Fx > 211 



would be then in agreement with the pulsar association of IFGL J1312. 6-1-0048 (see .Marelli et al.ll201 ih . 



173 3.2. IFGL J1739.4+8717 

174 We detected one X-ray source within the 2FGL error ellipse of IFGL J1739.4-I-8717 (designated as 

175 Suzaku J1742-I-8715). The detection significance is la (total of 204 net source counts from three detectors). 

176 The XIS image of Suzaku J1742-I-8715 seems to be relatively diffuse. However, Swift XRT has detected 

177 recently the same X-ray emitter as a point source at [RA, Dec] = [265.430, 87.245] with 90% position error 

178 radius rgo% = 6".9. Hence we conclude that in a relatively short Suzaku exposure (16.7 ksec) the apparently 

179 diffuse structure of the object is just an artifact of a low photon statistics. We found several weak optical 

180 sources coinciding positionally with the X-ray source, but no radio counterpart. 



^See Swift Technical Handbook (http : / /heasarc.nasa. gov/ docs /swift /proposals/ appendixj: . html) 
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181 We extracted the source photons from a 2' radius circle around Suzaku J1742+8715 and the background 

182 photons from a 3' radius circle. The resulting X-ray spectrum and X-ray light curve with the time bin of 

183 5760 s are shown in Figures|5] and |6l respectively. A constant fit to the light curve of the X-ray source 

184 returned x^/d.o.f. = 1.9/6, indicating that the flux was steady during the Suzaku exposure with the 

185 probability of > 93%. The spectrum was initially fitted by an absorbed power-law model (x^/d.o.f. = 

186 11.4/12) and this model returned Nn = 0. We then refit the spectrum after fixing the Nu value to zero and 

187 obtained a photon index of F = 2.1ibj;Q | and the unabsorbed 2 — lOkeV flux 3.6ti'5 cm^^ s^^. 

188 We found one relatively bright radio source inside the 2FGL error region of IFGL J1739.4-I-8717, 

189 namely NVSS J173722-I-871744 (~ 60 mJy at 1.4 GHz). Its position is marked with a g reen rhombus in Fig- 



191 



200 



urelH This radio sour ce is detected at several frequencies (151 MHz and 325 MHz; see lfialdwin et al.lll985 



Rengehnk et al.lll997L respectively), but has no obvious X-ray counterpart in our Suzaku data. Assuming a 

192 power-law form of the radio continuum (Fjy oc z^~"0> we calculate the spectral index of NVSS J173722-I-871744 

193 as ttr = 0.2 ifc 0.1 from the archiv al data. This is, in fact, a typical radio spectral index for blazar sources 



194 (e.g.. lSowards-Emmerd et al.ll2005h . Moreover, we also found an optical and infrared counterpart, USNOB 

195 1772-0020476 ([RA, Dec] = [264.35499(2), 87.29532(2)]), which is located only 2".2 away from the NVSS 

196 source with optical magnitudes, B2 = 19.30, R2 = 17.61, and I = 17.07, and 2MASS J17372480-I-8717433 

197 ([RA, Dec] = [264.3533(1), 87.2953(1)], and magnitudes J = 16.2, H = 15.5, and K = 15.8). This allowed 

198 US to calculate radio-to-optical and optical-to-X-ray spectral indices as Oro = 0.40 ± 0.01 and aox > 1-77, 

199 with the latter utilizing the X-ray flux upper limit (see Figure [T6l). The se indices are cons istent with those of 



so-called intermediate synchrotron peaked blazars (see Fig. 7 in lAckermann et al. Il201lb . 



201 3.3. IFGL J2043.2+1709 

202 We found one X-ray point source within the 2FGL error ellipse of IFGL J2043.2-I- 1709 at the respective 

203 position [RA, Dec]=[310°. 801(2), 17°. 171(2)], designated as Suzaku J2043-I-1710, and a second source lo- 

204 cated on the edge of the ellipse at [RA, Dec]=[310°. 822(3), 17°. 133(3)], designated as Suzaku J2043-I-1707 

205 (hereafter srcA and srcB for short, respectively). These two X-ray objects are not recorded in any available 

206 X-ray source catalogs. Swift X RT has also detected the same two X-ray emitters. Moreover, Swift Ultravi- 



207 olet/Optical Telescope (UVOT; iRoming et al. Il2005h discovered a possible UV counterpart of srcA at [RA, 



208 Dec] = [310.80314(6), 17.17286(7)], for which we found an association with the optical source USNOB 

209 1071-0645302 ([RA, Dec] = [310.80315(7), 17.17300(4)], and magnitudes B2 = 19.89, R2 = 20.6). No ra- 

210 dio and infrared counterparts for srcA or srcB were found. The X-ray image of the targeted field is shown in 

211 FigureEl The radio position of the MSP PSR J2043-I-171 1 recently associated with the IFGL J2043. 2-1-1 709 

212 is marked in the figure with a green cross centered at [RA,Dec]=[3 10°. 8370129(2), 17°. 1913744(3)]. 

213 To extract the spectra and light curves of the detected X-ray sources, we set the source regions and the 

214 background region as indicated in Figure|7] The detection significances of both sources were calculated to 

215 be 16(7 (total of 397 net source counts from three detectors) and 8a (total of 155 net source counts from three 

216 detectors), respectively. No X-ray counterpart of PSR J2043H-1711 was detected in our Suzaku exposure. 
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217 
218 



223 



231 



The light curves of srcA and srcB with time bins of 5760 s are given in FigureHl and the corresponding 
spectra are presented in Figure|9] Both light curves can be well fit by constant count rates (x^/d.o.f. = 8.5/9 



219 and X /d.o.f. = 2.6/9 for srcA and srcB, respectively). 

220 An absorbed power-law model provided the best fit to the spectrum of srcA, returning the photon index 

221 r = 1.67+0^5 and x^/d.o.f = 20.4/19. In the fit, the value of Nn was frozen at the Galact ic value in the 



direction of the target taken from the LAB Survey of Galactic HI, i.e., 6.63 x 10^*^ cm ^ (IKalberla et al 
2005b . We note that the fit with the absorption column density set free returned the value for A^h consistent 

224 with the Galactic one. The derived unabsorbed X-ray flux of srcA in the 2 — lOkeV photon energy range 

225 is 2.1^04 X 10^^^ ergcm^^ s^^. Similarly, the X-ray spectrum of srcB was initially fitted with an absorbed 

226 power-law model. However, this model fit returned negligible value of Nh and therefore we fixed A^h = 0. 

227 The best fit model with x^/d.o.f. = 8.6/10 returned then the photon index F = 1 •61^0 3 and the unabsorbed 

228 2 - 10keVflux9.ll|^ x IQ-i^ erg cm^^ g-i. 

229 Assuming a power-law model with photon index F = 2, we calculated upper limits (90% confidence) 

230 in the 0.5 - 2keV and 2 - lOkeV bands at the position of PSR J2043-I- 171 1 as < 2.4 x 10"^'' erg cm^^ g-i 
and < 3.6 x 10^^^ ergcm^^ s^^, respectively. This gives the ratio of the 0.1 — 100 GeV an d 2 — lOkeV 



232 fluxes F^/Fx > 864 which is consistent with the values claimed for the GeV-detected MSP (IMarelli et al. 



201 ih . 



234 3.4. IFGL J2302.8+4443 

235 One X-ray point source was discovered inside the 2FGL error ellipse of IFGL J2302. 8-1-4443 at the 

236 position of [RA, Dec]=[345°. 695(3), 44°.707(2)] (designated as Suzaku J2302-I-4442). The detection signif- 

237 icance of this X-ray source was calculated as 5.35cr (total of 121 net source counts from three detectors). As 

238 shown in Figure[TOl the detected X-ray source coincides with the radio position of the MSP PSR J2302-I-4442 

239 (green cross centered at [RA, Dec]=[345°. 695748(3), 44°.706136(1)]). This MSP was recently claimed to 

240 be associated with IFGL J2302. 8-1-4443. No optical and infrared counterparts of the pulsar were found. In 

241 the analysis of the Suzaku data, the source and the background regions were set as indicated in the Figure. 

242 The light curve of the X-ray counterpart with the time bin of 11520 s is presented in Figure[TT] The applied 

243 test assuming a constant count rate gave x^/d.o.f= 4.98/5, indicating that the X-ray flux was steady during 

244 the Suzaku exposure. 

245 The XIS spectrum of the detected X-ray source is shown in Figure[T2l Initially, we fitted the spectrum 

246 with an absorbed power-law model. This model returned an extremely soft continuum characterized by a 

247 photon index of F = 4.5lf y, and a rather high value of A'^h = 0.56to 34 x 10^^ cm^^. This absorption 

248 value is in exc ess of the Galactic col umn density in the direction of the target, namely 0.132 x 10^^ cm~^, as 



determined by IKalberla et al.l(|2005h from the LAB survey of Galactic HI. Next we applied an absorbed black 



body model but the model fit returned Nu consistent with zero. We therefore fixed A'h = and obtained the 
best-fit (x^/d.o.f. = 3.0/11) temperature of kT ~ O.Slto ogkeV (see Figure[l2ll. The observed X-ray flux 
of the source is 2.61q5 ^ 10^^'^ erg cm^^ s^^ in the 0.5 — 2keV range and 3.3 ^20 ^ 10^^^ erg cm^^ s^^ 
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253 in the 2 - lOkeV range. The ratio between 7-ray and 2 - lOkeV fluxes for IFGL J2302.8+4443/Suzaku 

254 J2302+4442 reads as Fj/Fx = 14531. This ratio is relatively high but still consistent with those of Fermi- 

255 LAT /LAT pulsars reported in lMarelli et all (120111) . 



256 



3.5. Other Sources 



257 
258 
259 
260 
261 
262 
263 
264 
265 
266 
267 
268 
269 
270 
271 
272 
273 
274 



No X-ray sources were detected in our Suzaku exposures inside the 2FGL error ellipses of the re- 
maining 7-ray objects: IFGL J0106.7-F4853, IFGL J1743. 8-7620, and IFGL J1902.0-5110. The cor- 
responding X-ray images of the targets are shown in Figures[T3l [T4j and [15] Out of these three sources, 
IFGL J1902.0-5110 has been recently associated with the MSP PSR J1902-5105. Even more recently, 
after the submissi on of this paper, also IFGL J0106.7-I-4853 has been identified with a new 7-ray pulsar 
PSR J0106-I-4855 (jpietsch et al. Iboilk Thus only IFGL J1743. 8-7620 remains at the moment unasso- 
ciated. We calculated an X-ray upper limit in the 2 — 10 keV band for the latter source assuming X-ray 
emission from a point source located at the center of 2FGL error region, obtaining in this way the energy 
flux < 4.5 X 10~^^ erg cm^^ s^^. In the case of IFGL J1902.0— 5110, we set the X-ray extraction region as 
a 1' radius circle around the radio position of PSR J1902— 5105. The resulting X-ray upper limits are calcu- 
lated as < 2.0xl0"i'^ergcm-2s-i (0.5-2 keV) and < 2.5 x lO^^"^ erg cm"^ ^-i (2-10keV)by assuming 
an absorbed power-law model with photon index F = 2. Finally, X-ray upper limits for IFGL JO 106.7-1-4853 
were calculated for a circular extraction region with radius 1' around the position of PSR J0106-I-4855, ob- 
taining < 7.7 X 10~^^ergcm"2s-i (0.5 - 2keV) and < 1.0 x 10"^^ erg cm"^ ^-i (2 - 10 keV). For these 
two pulsars, the ratio of the 0.1 - 100 GeV and 2 - lOkeV fluxes are F^/Fx > 976 and F^ lFx > 2640, 
respec tively. These flux ratios are consistent with the ratios of Fermi-L AT pulsars reported in lMarelli et all 
(|201lh . The results of the Swift observations of all of the three targeted 7-ray emitters are consistent with 
the Suzaku results. 



4. Discussion and Conclusions 



276 
277 
278 
279 
280 
281 
282 
283 
284 



In this paper, we report on the results of X-ray follow-up observations of seven bright Fermi-hAI 
sources at high Galactic latitudes {\b\ > 10°) using Suzaku XIS. We discovered the X-ray counterpart 
of IFGL J2302.8-I-4443 coinciding with the position of the MSP PSR J2302-I-4442 recently claimed to 
be associated with the 7-ray emitter. We did not however, detect X-ray counterparts for the other four 
Fermi-LAI objects similarly identified with a normal pulsar (IFGL J0106.7-I-4853) and MSPs, namely for 
IFGL J1312.6-F0048, IFGL J1902.0-5110 and IFGL J2043.2-F1709. (In a few cases the X-rays sources 
have been detected within the 2FGL error ellipses, but none at the positions of the pulsars.) A relatively 
weak X-ray source was found inside the 2FGL error region of IFGL J1739.4-I-8717. Finally, no candidate 
for the X-ray counterpart was detected for the remaining object IFGL J1743.8— 7620. 



285 Including ourprevious observati ons of IFGL J1231. 1-1410, IFGL J131 1.7-3429, IFGL 11333.2-^5056, 

286 and IFGL J2017. 3-1-0603 reported in iMaeda et al.l (|201lh . our sample of high Galactic-latitude Fermi-LAI 
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287 objects initially selected as unidentified and studied with Suzaku consists now of eleven targets. For eight 

288 of these, we have detected single or multiple X-ray sources within the LAT error ellipses. Over the time 

289 period when the Suzaku observations were being obtained, six targets from the 7-ray sample were found 

290 to be associated with MSPs, one target (IFGL J0106.7+4853) has been associated with a normal pulsar, 

291 and one source (IFGL J1333. 2+5056) has been classified as an AGN candidate, all in agreement with the 

292 gathered X-ray data. Still, four objects from the list remain unidentified. As argued below, one of these four, 

293 IFGL J1739.4-I-8717, is quite likely a high-redshift blazar. 

294 The source IFGL J1739.4-I-8717 was characterized by an enhanced flux level within the LAT photon 

295 energy range during the the first seven months of Femii-LAI operatior|H After that time, the activity of the 7- 

296 ray emitter decreased. The p hoton index in the LAT energy band is = 2. 1 ± 0. 1 (where dN/dE cx E"'"^^ is 

297 the differen tial photon flux; ^Abdo e t aPboil ), which is a typical value for the 7-ray spectra of BL Lac type 

298 blazars (see lAckermann et al. .2011.') . Importantly, as written in section IJ!2] one relatively bright radio source, 

299 NVSSJ173722-I-871744, is located inside the 2FGL error region of IFGL J1739.4-F8717 (see FigureSll. 

300 With a typical radio spectral index for blazar sources and radio-to-optical and optical-to-X-ray spectral 

301 indices that are consistent with blazar broadband spectrum, it is quite likely that IFGL J1739.4-I-8717 is 

302 indeed associated with a distant blazar currently characterized by an activity level low enough so that its 

303 X-ray emission was below the detection limit of the XIS instrument (~ 10^^^ erg cm" 

304 the performed Suzaku observations. 



313 
314 
315 
316 
317 
318 
319 
320 
321 



^ s ^) at the time of 



In general, unidentified sources constituted a large fraction of the population of 7-ray emitters detected 
by EGRET (~ 60% in 3EG), and at present about 31 % of Fermi-LAI sources in 2FGL catalog remain 
unassociated (specifically 273 sources at high Galactic latitudes (|6| > 10°) and 303 sources at low Galactic 
latitudes \b\ < 10°). Those located at the lowest Galactic latitudes (\b\ < 5°) are most widely expected to be 
associated with local systems such as molecular clouds, sup ernova remnants, massive stars, high-mass X-ray 
binar i es, radio quiet pulsa rs, and pulsar wind nebulae (e.g.. lKaaret & Cottamlll996l : lYadigaroglu & Romani 
1997l : lRomero et al. II1999I) . I n particular, half a doze n of the brightest 3EG sources in the Galactic plane were 



identified as young pulsars ( Thompson et al. 19991 ). despite the relatively poor localization of the EGRET 
sources and the source confusion complicated substantially the identification procedure. On the other hand, 
most of the unassociated 3EG sourc es at high Galactic latitudes (|6| > 10°) were later identified as blazars 
(|Sowards-Emmerd et al l2004l ). Pulsars were therefore expected to be found mainly among GeV emit- 
ters at low Galactic latitudes, while blazars were supposed to constitute the main population of GeV emitters 
at high Galactic latitudes. But the 7-rav-bright pulsa rs were also expected to be found at intermediate Galac- 
tic latitudes (5° < |6| < 73° : ICrawford et al.ll2006h . The identification of a number of Fermi-LAI objects 
located above the Galactic plane (\b\ > 10°) with such systems (pred ominantly with MSPs) confirmed 



these expectations (see, e.g.. Ransom & Pulsar Search Consortium 201(1 and the discussion in lMaeda et al. 



2011 



322 In FigurefTT] we plot the X-ray-to-7-ray energy flux density ratios (-F2-iokcv/^o. 1-100 Gev) versus 

323 radio-to-7-ray energy flux density ratios (-Fi.4GHz/-^o.i-iooGcv) for theFerm/-LAT objects from our Suzaku 



http:// he asarc.gsfc.nasa. gov/FTP/ f ermi/ data/lat / catalogs / source /light curves /2FGLJ1 7 38 . 9 + 8716. png 
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338 
339 
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sample discussed h ere (blue circles) with the radio data available in the literature (lAbdo et all l2010bl : 



350 
351 

352 
353 
354 
355 
356 
357 
358 



Condon et alJll998L see Table|3]). These can be compared with the analogous ratios evaluated for bright 



Fermi-LAI objects identified with blazars and MSPs (denoted in the figure by red crosses and squares, re- 
spectively). We remind the reader that the blazar class includes flat spectrum radio quasars (FSRQs) and 



BL La certae objects (BL Lacs). In addition, in the figure we plot the two targets discussed in iMaeda et al. 



(|201lh . namely IFGL J1333.2+5056 most likely associated with an AGN, and apeculiar object IFGL J1311.7-3429 
(pink stars). As shown, the blazar and MSP populations are clearly separated in the constructed flux ra- 
tio plane. Also, four objects from our sample which have recently been associated with a normal pulsar 
(IFGL J0106.7-F4853) and MSPs (IFGL J13 12.6-^0048, IFGL J2043.2-F1709, and IFGL J2302. 8-^4443) 
occupy the same region in the analyzed parameter space as the previously known MSPs detected in the GeV 
rangeQ On the other hand, in the case of IFGL J1739.4-I-8717 the evaluated energy flux density ratios — 
which are very similar to those characterizing IFGL J1333. 2-1-5056 — are consistent with the blazar iden- 
tification proposed above, if only NVSS J173722-I-871744 is considered as the true counterpart of the 7-ray 
emitter. 

In all, we conclude that the gathered Suzaku XIS data together with the broad-band properties of 
the analyzed Fermi-l^P^ objects are in agreement with the identification of most of them as MSPs. Yet 
a few cases in the analyzed sample (IFGL J1739.4-I-8717, IFGL J1333. 2-1-5056) constitute quite proba- 
ble associations with AGN (high-redshift blazars). Finally, the nature of the remaining few targets (like 
IFGL J1311.7— 3429) is still an open question, although, as inferred from Figure [TtJ the MSP identifica- 
tion seems more viable than the blazar one. In the near future, we are further continuing our X-ray studies 
during the Suzaku A06 cycle, focusing on both new targets (IFGL J0I03.1-i4840, IFGL J1946.7-5404, 
and IFGL J2339.7— 0531), but also performing ultra-deep exposures on particularly intriguing sources like 
IFGL J1311.7-3429. 

Finally, let us comment in more detail on the case of IFGL J2302. 8-1-4443. This object, as already 
mentioned above, has recen tly been associated w ith the millisecond pulsar PSR J2302-I-4442 discovered by 
the Nan9ay radio telescope ( Cognard et ahlboil ). The rotation period of the pulsar is P ~ 5.19 ms, the spin- 
down luminos i ty is E ~ 3.74 x lO^^ergs"^, and the characteristic age can be evaluated as r ~ 6.2 Gyr. 



Cognard et al.l (|2011h reported also on the detection of the X-ray counterpart of the pulsar with XMM- 

-0.4 

1.1 



Newton, with the unabsorbed 0.5 — 3keV flux of 3.11^4 x 10 erg cm ^ s ^ . This is consistent with our 



Suzaku detection (the re-calculated flux in the same photon energy range 2.9'^\^2 ^ 10^^^ erg cm~^ s~^). 
Anticipating the dista nce of the pulsar d 1. 18 kpc which was inferred from the NE2001 model of Galactic 
free electron density dCordes & Lazid 120011) . the X-ray luminosity of PSR J2302-I-4442 can therefore be 
evaluated roughly as Lx ~ 3 x 10^" erg s~^. This, together with E provided above, is then in good agreement 
with t he scaling relation between the X-ray and spin-do\yn luminosities ~ 10^^ x E established for 
MSPs feecker & Truempeilll997l : lGaensler & SlanelboodJzhang et aLll2007h . although the inferred distance 



"'in the case of IFGL J1312.6+0048, the X-ray-to-7-ray energy flux density ratio shown in FigurefTTlis evaluated assuming the 
association of the 7-ray sources with the MSP PSR J1312-I-00. That is, the Suzaku XIS upper limit derived at the position of the 
pulsar is considered, and not the X-ray flux of the Suzaku source detected within 2FGL error region. 



-12- 



359 
360 

361 
362 
363 
364 
365 
366 
367 
368 
369 
370 

371 
372 
373 



of PSR J2302+4442 is indicated to be smaller by a factor of four consideri ng an unphysically h igh 7-ray 
efficiency and, instead, assuming average efficiency of 7-ray MSPs ~ 10 % (iCognard et alj|201ll) . 

The spectrum of the X-ray counterpart of PSR J2302-I-4442/1FGL J2302. 8-1-4443 was well fit by a 
blackbody model, and this is again in agreement with the idea that the observed X-ray photons originate 
from thermal emission from the surface of a rotating magnetized neutron star. Because of the limited photon 
statistics in the higher energy range of XIS, we could not however confirm the presence of a non- thermal 
component above 2keV (clearly detected in the case of IFGL J1231. 1 — 1410 by iMaeda et al.ll201 ih . Inter- 
estingly, for the isolated pulsars as old as PSR J2302-I-4442 (c haracteristic age of about 6 Gyr), the surface 
tempe rature of the neutron star is expected to be T < 10^ K (INomoto & Tsurutalll987l: IPage & Applegate 
1992h . Our detection of the X-ray counterpart indicated /cT ~ 0.31 ± 0.03 keV, i.e. T ~ 3.6 x 10^ K instead. 



Some reheating pr ocess is there fore required, possibly related to the impact of relativistic particles on polar 
caps ( Becker & Truempenl 19971. and references therein). 

C.C. Cheung's work at NRL is sponsored by NASA DPR S-15633-Y. L.S. is grateful for the support 
from Polish MNiSW through the grant N-N203-380336. We thank the anonymous referee for a careful 
reading of the manuscript and useful suggestion which helped to improve the paper. 
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Table 1 : Suzaku Observation Logs and 7-ray properties of the targets 



Name 


DBS ID 


Pointing Direction" 


Observation 


Effective 


7-ray photon 


7-ray 


flux'' 










start 


exposure 


index'' 


-fo.l-lOOGcV 






RA [deg] 


DEC [deg] 


(UT) 


[ksec] 




[erg cm 




IFGL J0106.7+4853 


705010010 


16.6433 


48.9425 


2010/07/15 09:53:44 


19.7 


2.05 


2.64 X 


10-" 


IFGL J1312.6+0048 


705011010 


198.1860 


0.8370 


2011/01/17 04:38:43 


17.5 


1.99 


1.96 X 


10-" 


IFGL J1739.4+8717 


705012010 


264.8730 


87.2900 


2010/04/26 23:41:56 


16.7 


2.12 


2.92 X 


10-" 


IFGL J1743. 8-7620 


705013010 


265.9610 


-76.3420 


2010/04/14 00:16:05 


34.0 


2.14 


2.83 X 


10-" 


IFGL J1902.0-5110 


705014010 


285.5220 


-51.1700 


2010/04/13 06:37:43 


38.5 


2.10 


2.46 X 


10-" 


IFGL J2043.2+1709 


705015010 


310.8220 


17.1640 


2010/05/03 19:23:03 


17.7 


2.13 


3.15 X 


10-" 


IFGL J2302.8+4443 


705016010 


345.7070 


44.7230 


2010/06/26 05:52:36 


35.6 


2.04 


4.81 X 


10-" 



" The pointing directions are the values of pl anned target coordi nates. 
These values are taken from IFGL catalog jAbdo et alJboiObh . 



Table 2: Results of the Suzaku observations of the selected Fermi-L AT sources. Positions of detected X-ray 
sources are presented. 



Name 


Position 




Suzaku detection 


RA[deg] 


DEC [deg] 


IFGL J0106.7-^4853 








IFGL J1312.6-^0048 


SUZAKU 11312-^0050 


198.235(2) 


0.835(2) 


IFGL J1739.4-^8717 


source" 






IFGL J 1743. 8 -7620 








IFGL J1902.0-5110 








IFGL J2043.2-I-1709 


SUZAKU 12043-^1710 


310.801(2) 


17.171(2) 




SUZAKU 12043-^1707 


310.822(3) 


17.133(3) 


IFGL 12302.8-1-4443 


PSR J2302-I-4442 


345.695(3) 


44.707(2) 



"We could not determine the exact position of the source because of extended source image (see section 3.2). 
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Table 3: Results of the analysis of the Suzaku fXlS data and radio fluxes of the sources 



Name 


Nh 


Model 


Model 


X-ray flux 


Radio flux 










parameter 


F2-10keV 


Fl.4GHz 






[10^° cm"^] 






[ergcm"^ s~^] 


[mJy] 


IFGL J0106.7+4853 


PSR JO 106+4855 


11.9 (fixed) 


PL 


r=2.0 (fixed) 


< 1.0 X 10~" 


< 0.008 


IFGL J1312.6+0048 


source 


0.0 (fixed) 


PL 




S.Otli X 10-" 


<0.71 




off-source'' 


2.1 (fixed) 


PL 


r=2.0 (fixed) 


< 9.3 X 10"" 


<0.82 


1FGLJ1739.4+8717 


source 


6.36 (fixed) 


PL 


p_7 1 +0.5 
-^■'--0.4 


3.6 tli X 10-" 


<0.58 




NVSS J 173722+87 1744 


6.4 (fixed) 


PL 


r=2.0 (fixed) 


< 1.5 X 10"" 


61.3 


IFGL J1743.8-7620 




8.13 (fixed) 


PL 
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IFGL J 1902.0-5 110 


PSR J 1902-5 105 


4.87 (fixed) 


PL 
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src A 


6.63 (fixed) 


PL 


P_i 67+0.16 
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< 1.01 




src B 


0.0 (fixed) 


PL 


r-1 6+°* 

1 -i.O_g 3 
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< 1.00 
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6.67 (fixed) 


PL 


r=2.0 (fixed) 


< 3.6 X 10"" 


0.01" 


IFGL J2302.8+4443 


PSR J2302+4442 


0.0 (fixed) 


BB 


fcT=0.3lin5keV 


3.3l2:o X 10"" 


1.2'' 



Note. — Summary of analysis results of the seven studied IFGL sources. Unabsorbed X-ray flux and 90 % confidence upper limits 
in the 2—10 keV band for all the sources are listed. The X-ray upper limits are calculated assuming a power law model with a 
photon index of 2 for the IFGL sources with no X-ray detection. In the seventh column, radio fluxes or upper limits at 1.4 GHz 
are presented. Radio fluxes and upper limits The radio values with no footnote marks were taken either directly from the NVSS 
(Condon et al. 1998) catalog (for NVSS J173722+871744) or measured by us from the NVSS images (the upper limits; 90% 
confidence levels). 

" This value is estimated by averaging the radio flux over four observations during which the pulsar was detected. 
From Cognard et al. (201 1). 

The evaluated off-source X-ray upper limits expected to represent MSP PSR J1312+00. 
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Fig. 1.— X-ray image of IFGL J1312.6+0048 by SuzakufXlSO+3 (FI CCDs) in the 0.5 to 10 keV energy 
band. Thick solid ellipse denotes the 95% position error of IFGL J1312.6+0048 in the 2FGL catalog. Thin 
solid and dashed circles show the source and background regions, respectively. The accurate position of 
PSR J1312+00 is still not available in the hterature. 
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Fig. 2. — Suzaku/XlS light curve of the possible X-ray counterpart of IFGL J1312.6+0048. 
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Fig. 3. — Suzaku/XlS spectrum of the possible X-ray counterpart for IFGL Jl 3 12.6+0048 fitted with a 
power-law model. Black plots show the FI data and gray plots show the BI data. 
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Fig. 4.— X-ray image of IFGL J1739.4+8717 by SuzakufXlSO+3 (FI CCDs) in the 0.5 to 10 keV 
energy band. Thick solid ellipse denotes the 95% position error of the 2FGL catalog counterpart of 
IFGL J1739.4+8717. Thin solid and dashed circles show the source and background regions, respectively. 
The radio position of NVSS J173722+871744 is marked with a green rhombus. 
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Fig. 5. — X-ray spectrum of the X-ray counterpart of IFGL J1739.4-I-8717 fitted with a power-law model. 
Black plots show the FI data and gray plots show the BI data. 
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Fig. 6.— Suzaku /XIS light curve of the X-ray counterpart of IFGL J1739.4+8717. 



-23- 



17.400 



17.350 



17.300 



fl) 17.250 

0) 

2, 

S 17.200 

o 

o 

CM 



E 17.150 
o 



Si 17.100 
Q 



17.050 



17.000 




311.000 



310.900 



310.800 



310.700 



310.600 



Right Ascension (J2000) [degree] 



OOe+OO 2.99e-08 5.98e-08 9.00e-08 1.20e-07 1.50e-07 1.80e-07 2.10e-07 2.40e-07 2.70e-07 3.00e-07 



Fig. 7.— X-ray image of IFGL J2043.2+1709 by SuzakufXlSO+3 (FI CCDs) in the 0.5 to 10 keV energy 
band. The thick sohd ellipse denotes the 95% position error of IFGL J2043.2+1709 from the 2FGL catalog. 
Thin solid and dashed circles show the source and background regions, respectively. Two X-ray sources 
inside the positional error ellipse of IFGL J2043. 2-1-1 709 are named srcA for the northern source and srcB 
for the southern source. The position of the associated MSP PSR J2043-I-171 1 is marked with a cross. 
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(b) Source B 

Fig. 8. — Suzaku/XlS light curve of srcA (a) and srcB (b). 
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Fig. 9. — SuzakufXlS spectra of srcA (a) and srcB (b) inside the error ellipse of IFGL J2043.2+1709. Both 
spectra are fitted with an absorbed power-law model. Black plots show the FI data and gray plots show the 
BI data. 
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Fig. 10.— X-ray image of IFGL J2302. 8+4443 by SuzakufXlSO+3 (FI CCDs) in the 0.5 to 2 keV en- 
ergy band. Thick solid ellipse denotes 95% position error of IFGL J2302. 8-1-4443 from the 2FGL catalog. 
Thin solid and dashed circles show the source and background regions, respectively. The position of the 
associated MSP PSR J2302-H4442 is marked with a cross. 
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Fig. 11.— SuzakufXlS light curve of the X-ray counterpart of IFGL J2302.8+4443. 
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Fig. 12. — Suzaku/XlS spectrum of the X-ray counterpart of IFGL J2302.8+4443 fitted with a black body 
model. Black plots show the FX data and gray plots show the BI data. 
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Fig. 13.— X-ray image of IFGL J0106.7+4853 by SuzakufXlSO+3 (FI CCDs) in the 0.5 to 10 keV energy 
band. Thick solid ellipse denotes 95% position error of IFGL J0106.7+4853 from the 2FGL catalog. 
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Fig. 14.— X-ray image of IFGL J1743.8-7620 by SuzakufXlSO+3 (FI CCDs) in the 0.5 to 10 keV energy 
band. Thick sohd ellipse denotes 95% position error of IFGL J1743.8— 7620 from the 2FGL catalog. 
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Fig. 15.— X-ray image of IFGL J1902.0-5110 by Suzaku/XlSO+3 (FI CCDs) in the 0.5 to 10 keV energy 
band. Thick solid ellipse denotes 95% position error of IFGL J1902.0— 5110 from the 2FGL catalog. The 
accurate position of PSR J1902— 5 105 is still not available in the literature. 
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Fig. 16. — Spectral energy distributions of IFGL J1739.4+8717. The X-ray data represents 90% confi- 
dence upper limit c alculated from our Suzaku data. The 7-ray data points are taken f rom the IFG L catalog 
jAbdo et al.lboiObI '). The radio data points a re 6C B 175708+8719 24 in 6C catalog fealdwin et al. 1985), 
WENSS B1758.4+8718 in W ENSS catalog jRengelink et al.ll 19971^ and NVSS J173722+871744 in NVSS 
catalog dCondon et al.lll998h . Infrared and optical plots, 2 MASS J17372480+87 17433 and USNOB 1772- 
0020476, are quote d from 2MASS point source catalog dskrutskie et al. 1120061) and USNO Bl.O catalog 
dMonet et al. 112003b . respectively. 
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Fig. 17. — X-ray to 7-ray flux ratios versus radio to 7-ray flux ratios for unassociated Fermi-hAI sources. 
The blue circles show the flux ratios of unassociated Fermi-hAI sources which we observed using the 
Suzaku satellite. The 7-ray flux and the radio upper limits of these sources are taken from values of the 
IFGL catalog and the NVSS catalog (Condon et al. 1998), respectively. Red x marks indicate the data plots 
of blazars listed in Abdo et al. (201 Of). Red squares represent the data points of the Fermi-LAT sources that 
are associated with MSPs discovered by Cognard et al. (2011) and Ransom et al. (201 1), and three Ferrni- 
LAT sources listed in Abdo et al. (2010j) for which the best fit model of the X-ray spectra were determined. 
The values of X-ray flux and radio fluxes of MSPs in Abdo et al. (2010j) are quoted and derived from Bailes 
et al. (1997), Lundgen et al. (1995), Navarro et al. (1995), Sakurai et al. (2001), Webb et al. (2004a), and 
Webb et al. (2004b). For the MSPs we cannot quote the radio flux at 1.4 GHz, we estimated 1.4 GHz flux 
by the extrapolation using radio fluxes that are available in the literature assuming a power-law model. Pink 
stars show the two still unassociated IFGL sources from our first year campaign targets (Maeda et al. 201 1) 
and another two IFGL sources studied in Maeda et al. (2011) and that are now associated with MSPs are 
involved in red squares. 



